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^ . Abstract 

o . 

I We present the next-to-next-to-next-to-leading O(a^mq) result for the ground 

' state energy of a heavy quarkonium system. On the basis of this result we deter- 

■ mine the bottom quark mass from T(15) resonance and provide an explicit formula 

I relating the top quark mass to the resonance energy in ti threshold production. 
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o : 1 Introduction 

' Quark masses enter the Lagrangian of Quantum Chromo dynamics (QCD) as fundamental 
^ ■ parameters. Thus, it is of primary interest to perform accurate determinations of the mass 
values combining experimental input with precision calculations from the theoretical side. 
In this letter we focus on the bottom and top quark masses. The bottom quark mass 
is of particular interest due to the current and future B physics experiments. In the 
observables used to extract the Cabibbo-Kobajashi-Maskawa (CKM) matrix elements 
and to gain deeper insight in the understanding of CP violation the bottom quark mass 
enters as a crucial input quantity. Thus a precise knowledge of the latter is essential for 
the interpretation of the experimental data. On the other hand, the top quark mass is 
one of the key parameters for precision tests of the standard model of the electroweak 
interactions at high energies and for the search of "new physics" at a future e~^e~ linear 
collider. 

Whereas the masses of the up, down and strange quark require major non-perturbative 
input the higher quark masses can be computed with relatively high accuracy using per- 
turbative methods. The study of nonrelativistic heavy- quark- ant iquark systems is one 
of the most promising approaches in this respect. These systems allow for a model- 
independent perturbative treatment which relies entirely on first principles of QCD |jl|]. 
Nonperturbative effects are well under control for the top-antitop system and, at 



least for the ground state resonance, also for bottomonium. This makes heavy-quark- 
antiquark systems an ideal laboratory to determine the heavy-quark masses. Besides its 
phenomenological importance, the heavy-quarkonium system is very interesting from the 
theoretical point of view because it possesses a highly sophisticated multiscale dynamics 
and its study demands the full power of the effective-theory approach. 

Among other characteristics of the heavy quarkonium system the ground state energy 
is of primary phenomenological interest. It can be directly used to extract the bottom 
quark mass from the lowest resonance of the T family and essentially determines the res- 
onance energy in ti threshold production. The (^(afmg) heavy quarkonium spectrum was 
obtained first in 0] and afterwards confirmed in At (^(a^mg) the heavy quarkonium 

spectrum was derived in for vanishing /3-function and the large-/?o result can be found 
in Bi. 

In the present letter we complete the result of [0 for the ground state energy by 
calculating the remaining terms proportional to /3-function. We apply the result to the 
determination of the bottom quark mass from the mass of T(IS') resonance and give 
an explicit formula relating the top quark mass to the resonance energy in ti threshold 
production both for e~^e~ annihilation and 77 collisions. The main results are given by 
Eqs. (H), (H) and (|lD. 



2 Perturbative heavy quarkonium ground state en- 
ergy at 0{almq) 

For the calculation of the ground state energy we follow the general approach of p|JTO|. 
It is based on the nonrelativistic effective theory concept ||ll| , p!2| in its potential NRQCD 
(pNRQCD) incarnation ||l3l and is implemented with the threshold expansion [l^. In 



pNRQCD the dynamics of the nonrelativistic potential heavy quark- ant iquark pair is 
governed by the effective Schrodinger equation and by its multipole interaction to the ul- 
trasoft gluons. The relativistic effects of the harder modes that have been "integrated out" 
are encoded in the higher- dimensional operators which appear in the effective Hamiltonian 
and correspond to an expansion in the heavy quark velocity v, and the Wilson coefficients 
which are series in the strong coupling constant a^. 

The analysis of the heavy quarkonium spectrum at (9(afmq), i.e. the third order 
corrections to the Coulomb approximation involves two basic ingredients: the effective 
Hamiltonian and the retardation effect associated with the emission and absorption of 
dynamical ultrasoft gluons. The general form of the Hamiltonian valid up to next-to- 
next-to-next-to-leading orderQ (N^LO) is given in 0. A detailed discussion of the effects 
resulting from the chromoelectric dipole interaction of the heavy quark-antiquark pair to 



^In such a framework one has two expansion parameters, and v. The corrections are classified 
according to the total power of ag and v, i.e. N'^LO corrections contain terms of O^a^v^), with 
I + m = k. This has the consequence that, in general, different loop orders, which are counted in 
powers of as, contribute to the N'^LO result. Note, that for the corrections to the parameters of a heavy 
quarkonium bound state, where v ^ as, the standard ordering in powers of as is restored. 
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the ultrasoft gluons relevant for perturbative bound-state calculations at N^LO can be 
found in [[l5| ,[f|. 

For vanishing angular momentum we can write the energy level of the principal quan- 
tum number n as 

= E^n + 5EP + 5E^^ + 5E^^ + • • • , (1) 
where the leading order Coulomb energy is given by 
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SE^"* stands for corrections of order aj^. The first and second order corrections can be 
found in [^^0 for arbitrary n. For the ground state one gets 
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where = (Cpdsfng)), C(3) = 1.202057. . . is Riemann's (^-function, Cp = 4/3 and 
Ca = 3 are the eigenvalues of the quadratic Casimir operators of the fundamental and 
adjoint representations of the colour gauge group, respectively. The modified minimal- 
subtraction (MS) scheme for the renormalization of as is implied. The coefficients Oj 
parameterize the z-loop perturbative corrections to the Coulomb potential. The analytical 
results for the one- [|l6|,|l7l and two-loop |[T8|JT9|JT0[| coefficients are 
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where Ti? = 1/2 is the index of the fundamental representation and ni is the number of 
light-quark flavors. For convenience of the reader we also provide the flrst three coefficients 
of the QCD /?-f unction 
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The O (a^) corrections to the energy levels arise from several sources: 

(i) matrix elements of the N^LO operators of the effective Hamiltonian between Cou- 
lomb wave functions; 

(ii) higher iterations of the NLO and NNLO operators of the effective Hamiltonian in 
time-independent perturbation theory; 

(iii) matrix elements of the N^LO instantaneous operators generated by the emission 
and absorption of ultrasoft gluons; and 

(iv) the retarded ultrasoft contribution. 

The contributions to E^^ for vanishing /9-function were computed for general in . 
The corresponding correction to the ground state energy gets contributions from all four 
sources and reads 
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where S is the spin quantum number and La^ = — InlCpas)- The terms proportional to 
ai correspond to the iterations of lower-order operators as mentioned in point (ii) of the 
above list. The last term in Eq. (P) corresponds to retarded ultrasoft contribution (cf. 
point (iv)). The corresponding QCD Bethe logarithms can be found in |I5|,|7| with the 
numerical result Lf = —81.5379. The logarithmic \ia{as) part of Eq. (P) is known from 
previous analyses |20,21|. At present, only Pade estimates of the three-loop MS coefficient 
03 are available 



For the bottom and top quark case they read 
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Note that at N^LO the Wilson coefficients in the effective Hamiltonian are infrared (IR) 
divergent. The corresponding divergence is canceled against the ultraviolet (UV) one of 
the ultrasoft contribution. In our analysis based on the threshold expansion dimensional 
regularization is used to handle the divergences. For convenience, the IR singular pole 
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part of the Wilson coefficients is subtracted according to the MS prescription, so that the 
coefficients, in particular a^, are defined in the MS subtraction scheme both for UV and 
IR divergences. For consistency, the UV pole of the ultrasoft contribution is subtracted 
in the same way 0. 

The third order corrections to the ground state energy proportional to the /3-function 
only get contributions from (i) and (ii). The computation of the contribution from (i) 
is rather straightforward. It includes the corrections due to the three-loop running of 
the Coulomb potential which we mark by C.r., and the one-loop running of the l/rriq- 
suppressed and the Breit potential (B.r.) of the following form 
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The contribution of the type (ii) can be computed using the general approach of 
(see also |^). It includes the corrections due to the iterations of the one- and two- loop 
running of the Coulomb potential which we mark by C.i., and the iteration of the kinetic 
energy correction, the 1/mg and the Breit potential with the one-loop running of the 
Coulomb potential {B.i.). The corresponding analytical expressions read 
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where C(5) = 1.036928 .... As a cross-check to our analytical calculation we computed the 
non-trivial ingredients by solving the eigenvalue problem numerically in the limit — > 
with the help of the program schroedinger .nb ^4 . 
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By adding (||), ( ]T0| ) and ( ]Tl| ) we obtain the third order correction to the ground 
state energy proportional to the /5-function 
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The terms proportional to /3q in Eqs. (|T2|) can be found in . 

The total result for the third order correction to the ground state energy is given by 
the sum of Eqs. (0) and ([T2|) 
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Adopting the choice /i^ = C pds^^i 3)1^ q one obtains for the bottom and top system (for 
S=l) in numerical form 
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where we have separated the contributions arising from 03 and (3^. The only unknown 
ingredient in our result for 6Ei^^ is the three-loop MS coefficient 03 of the corrections 
to the static potential entering Eq. (|^). Up to now there are only estimates based on 
Fade approximation which we will use in our analysis. However, we will show that our 
final result only changes marginally even for a rather large deviations of 03 from its Fade 
estimate. 



3 Bottom and top quark mass determination 

In this section we focus on the bottom and top quark mass determination. We fix the value 
of the strong coupling constant which enters our theoretical expression to af\Mz) = 
0.1185 ±0.002 |5|. 
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Bottom quark. The starting point for the determination of the bottom quark mass 
is its relation to the mass of the T(IS') resonance 



Mx(i5) = 2m, + E^- + , (15) 



with Mx(i5) = 9.46030(26) GeV ||2^. Here is the nonperturbative correction to 

the ground state energy. The dominant nonperturbative correction is associated with the 
gluon condensate contribution given by 0,0] 

5--^-E, = / = 60 ± 30 MeV . (16) 

1275 (Cpa^mJ^ ^ ^ 

For the numerical estimate we use a^"^-* (/is) = 0.304 corresponding to the soft normalization 
scale /is = CFOis{[J^s)'mh characteristic to the Coulomb problem and nih = 4.83 GeV which 
is obtained from Eq. (pTSj) neglecting the corrections to the Coulomb binding energy. 
This is consistent with neglecting the perturbative correction to Eq. (|16|). For the gluon 
condensate we adopt the value from {asC^'^G^^^) = 0.04(2) GeV^ The effect of 



higher dimension condensates on the ground state energy is small Note, however 



that the contribution of the gluon condensate to the energy level of principal quantum 
number n grows as and becomes unacceptably large already for n = 2. Thus the 
excited levels are very sensitive to the nonperturbative long-wave vacuum fluctuations of 
the gluon field and cannot be used for reliable predictions. 

Combining Eq. (|I^) with the result of the previous section for the perturbative ground 
state energy up to 0{mqal) we obtain the bottom quark mass as a function the renor- 
malization scale of the strong coupling constant normalization, /i, which is plotted in 
Fig. 0(a). For the numerical evaluation we extract a^'^^nih) with rrih = 4.83 GeV from its 
value at Mz using four-loop /3-function accompanied with three-loop matching]^. a^^\mf,) 
is used as starting point in order to evaluate a^J^\^) at N'^LO with the help of the fc-loop 
P function. From Fig. ^(a) we see that the dependence on the renormalization scale be- 
comes very strong below /i ~ 2 GeV which indicates that the perturbative corrections 
are not under control. However, even above this scale the perturbative series for the pole 
mass shows no sign of convergence. This means that one can assign a numerical value to 
the pole mass only in a specified order of perturbation theory. 

On the contrary, it is widely believed that the MS mass mb(/i) at the scale /i = mb(/i) is 
a short-distance object which has much better perturbative properties. Thus, it seems to 
be reasonable to convert our result for the pole mass into mfc(mb). The relation between 
rub and rubijnb) is known up to three-loop approximation [p9| , p0| , pl[1 and shows sizable 
perturbative corrections. For this reason we suggest the following procedure to take into 
account these corrections in a most accurate way. The idea is that from the one-parametric 
family of m6(/i) we can choose a representative corresponding to some scale fi* in such a 
way that 

mbin*) = nib. (17) 
^We use the package RunDec pa to perform the running and matching of as 
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For a given fixed-order value of tlie pole mass Eq. ( p!7[ ) can be solved for /i*. In particular, 
for a pole mass to N'^LO we use the fc-loop relation between the MS and pole mass in 
Eq. (p!?!). Afterwards m^jirni,) can be computed from rnh{fi*) solving the RG equation. 
The advantage of this approach is obvious: we use the finite order relation between MS 
and pole mass at the scale where they are perturbatively close while the large difference 
between rnbiffib) and rrib is completely covered by the RG evolution which can be computed 
with very high accuracy as the corresponding anomalous dimension is known to four-loop 
approximation ||32,33]. The only restriction on the method could be connected to the 
value of n*. It should be large enough to allow for a reliable use of the RG equation. As 
we will see this condition is satisfied in practice. 

In Fig. |l](b) our result for rnb{jnb) is plotted at NLO, NNLO and N^LO as a function 
of the normalization scale which is used to obtain the pole mass nib (cf. Fig. |l](a)). It is 
remarkable that close to /z = 2.7 GeV, which is consistent with the physically motivated 
soft scale //^ ~ 2 GeV, both the second and the third order corrections vanish. This fact 
is a rather strong indication of the convergence of the series for mbirfib)- The numerical 
values for m^, rfibijnb) at yU = 2.7 GeV as well as the values of /x* and af'{^*) are listed 
in Tab. ^ for the different orders. 

An alternative approach to extract mb{jnb) was suggested in |34] within the so-called 
upsilon-expansion. The basic idea is to correlate the approximation for the energy spec- 
trum and the one of the quark mass relation by matching the leading powers of (3q. As 
a consequence the fc*^-order correction to the spectrum goes along with the [k + l)-loop 
mass relation. In this way the large corrections proportional to leading powers of (3q, 
which are associated to the IR renormalon contribution, are canceled in mbirfib). In order 
to perform the analysis at N^LO the four-loop relation between the MS and pole mass is 
necessary for which we adopt the large-/3o approximation result P^,P|. One can expect 
that it is close to the exact answer as the large-/3o approximation also works well at three- 
loop order |^9|,^. Adopting this prescription we obtain mbirfib) as shown in Fig. |l|(c). 
Note that within the upsilon-expansion the "scale of best convergence" is absent. Though 
the NNLO correction vanishes at /i ~ 2.7 GeV as it does within our yU* prescription the 
N^LO one at this point is sizable as can be seen from Tab. ^ This fact can be explained 
by the observation that the hypothesis of the renormalon contribution dominance in the 
perturbative series for the energy levels does not work at N^LO. To be precise, if one 
sets /i = /is in Eq. ([T2|) one finds that the renormalon induced term proportional to (3^ is 
indeed large, however, it covers less than a half of the total third order contribution to 
El which can be seen from Eq. (|14D. As a consequence, one would have to assign a larger 
uncertainty from the unknown higher order corrections to the N^LO value given in the 
last row of Tab. |l]. 

The uncertainty in the obtained value of rfibijfib) originates from several sources which 
are listed in Tab. ^ Beside the uncertainties in af\Mz) and the gluon condensate one 
has to take into account the errors in the Fade estimates of the coefficient 03, the charm 
quark mass effects and the higher order contributions. Let us consider the corresponding 
errors in more detail. In the above analysis we neglected the charm quark mass both 
in the correction to the spectrum and in MS-pole mass relation. As an estimate of the 
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Figure 1: (a) Bottom quark pole mass, nib, as a function of the renormalization scale /i, 
used for the numerical evaluation of Ei where the short-dashed, long-dashed and solid line 
corresponds to the NLO, NNLO and N^LO approximations, (b) MS bottom quark mass 
^bi^b), as a function of the renormalization scale /i, which is used for the extraction of 
the pole mass (cf. (a)). For the conversion from the pole to the MS mass the method 
described around Eq. (|T7D is used, (c) like (b), however, the upsilon-expansion is used. 
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order 


mfe (GeV) 


m;,(mb)W (GeV) 


fi* (GeV) 




mb{mby''P'^ (GeV) 


LO 


4.749 


4.749 






4.324 


NLO 


4.865 


4.348 


2.498 


0.2629 


4.223 


NNLO 


5.020 


4.338 


2.194 


0.2816 


4.225 


N^LO 


5.235 


4.346 


1.886 


0.3000 


4.274 



Table 1: Results for the pole mass and the MS mass as computed in the jj* scheme 
as described in the text, fj, = 2.7 GeV has been adopted. For completeness also fi* 
and a^'^\^*) are presented. In the last row the values for rfibijnb) according to upsilon- 
expansion are given. 



induced error we use the available results incorporating the charm quark mass effects 
in the spectrum [^,0 and the quark mass relation [^ J9[]. The value of higher order 



contributions is difficult to estimate. The nice convergence of the series for rfibijnb) in the 
first three orders implies their small size. Commonly the uncertainties from the higher 
orders are related to the uncertainty induced by a variation of the normalization scale /z. 
We estimate it by allowing a variation of /i by ±400 MeV, which corresponds to twice the 
difference between the scales at which the NNLO and N^LO corrections vanish, around 
our optimal value /i = 2.7 GeV. The uncertainty related to the error in the Fade estimate 
for 03 can be written as (a3"''~^''/98 — 1) x 10 MeV and will completely disappear once the 
exact computation becomes available. For the numerical estimate we allow for a ±100% 
deviation of 03 from the Fade estimate. 

Finally, our prediction for the bottom quark MS mass reads 

mb{mb) = 4.346 ± 0.070 GeV , (18) 

where the uncertainties are added up in quadrature. 

Now we are in the position to compare our result to the estimates obtained with other 
methods. There are different approaches to determine the bottom quark mass from the 
bb pair production close to threshold which are quite complementary as they are based on 
different experimental and theoretical input. Three basic techniques can be distinguished: 

(i) the direct analysis of the T resonances spectrum; 

(ii) low-order moments sum rules; 

(iii) the nonrelativistic T sum rules. 

Each approach has its own advantages and shortcomings. They are briefly characterized 
below and the results of the most recent theoretical analyses based on (i) - (iii) are listed 
in Tab. |[ 

The direct analysis of the spectrum of T resonances |I| used in this letter suffers from 
large nonperturbative corrections. The exited levels are too sensitive to the nonpertur- 
bative dynamics and only the ground state can be used for reliable predictions. On the 
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source 


Snibimb) (MeV) 


5Eres (MeV) 


asiMz) 


+13 
-8 


+149 
-156 




+66 
-47 


±41 


as 


±10 


±22 


(GG) 


±15 




m.c 


±10 








±10 


total 


±70 


±163 



Table 2: The sources of uncertainties and corresponding errors in the bottom quark 
MS mass and the resonance energy in tt threshold production. The intervals of the input 
parameters variation are specified in the text. The total uncertainty is obtained by adding 
up the individual ones in quadrature where in the case of asymmetric errors the maximum 
is chosen. 



other hand currently the most accurate perturbative description has been achieved within 
this framework. 

Study of the low-order moments of the spectral density within the sum rules approach 
relies on the global quark-hadron duality. Nonperturbative effects are negligible. The 
theoretical analysis consists in fully relativistic fixed-order calculations which have been 



performed up to NNLO [Q. The method allows for a direct determination of the MS 
mass. It is, however, sensitive to the behaviour of the spectral density above the threshold 
which is not well known experimentally. 

The nonrelativistic T sum rules are "between" the approaches (i) and (ii). They 
operate with the high (of the order l/a^) moments of the spectral density, which are 
saturated by the nonrelativistic near-threshold region. The experimental input is given 
by the parameters of the T resonances which are known with high accuracy. On the 
theoretical side the nonperturbative effects are well under control, however, the Coulomb 
effects should be properly taken into account. By now the analysis has been performed 
up to NNLO [^ ,|6l,^,Pl,P| . The extension to N^LO is a challenging theoretical problem. 
In particular, it requires the calculation of third order corrections to the heavy quarko- 
nium wave function at the origin. The first steps in this direction has been performed 
in 0,|T3|,^|2|,|31 



As we see from the Tab. ^ the results of all approaches are in reasonable agreement. 
Top quark. Let us now turn to the top quark case where the experimental data are 
not yet available. Our goal is to present a formula which can be directly used for the top 
quark mass determination from the characteristics of the cross section of tt production 
near threshold. One has to distinguish the production in e^e~ annihilation where the 
final state quark-antiquark pair is produced with S = 1 and in (unpolarized) 77 collisions 
where the dominant contribution is given by the S'-wave zero spin final state. 

The nonperturbative effects in the case of the top quark are negligible, however, the 
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ypfpypri pp 


mpf hon 


order 




PP98 H 
MY99 § 
BS99 |41 
HOO 
KSOl |38 

POl 44 
this work 


nonrelativistic T sum rules 

)) 

n 

5) 

low moments sum rules 
spectrum, T(IS') resonance 

77 


NNLO 
N^LO 


4.21 ±0.11 
4.20 ±0.10 
4.25 ±0.08 
4.17 ±0.05 

4.191 ±0.051 
4.210 ±0.090 ±0.025 

4.346 ± 0.070 



Table 3: Recent results for the MS mass obtained in different approaches and approxi- 
mations. 



effect of the top quark width has to be taken into account properly [45] as the relatively 



large width smears out the Coulomb-like resonances below the threshold. The NNLO 



analysis of the cross section shows that only the ground state pole results in the 
well-pronounced resonance. The higher poles and continuum, however, affect the position 
of the resonance peak and move it to higher energy. As a consequence, the resonance 
energy can be written in the form 

Eres = ± Ef-'- + ^r^Eres • (19) 

To compute the shift of the peak position due to the nonzero top quark width we 



use the result of ^3] for the e~^e tt and 77 — tt cross sections at NNLO. Numerically 
we get for both processes 

5^'E,es = 100± 10 MeV. (20) 

This value is rather stable with respect to the variation of all input parameters of our 
analysis such as the top quark mass and width, the value of the strong coupling constant 
and its renormalization scale. In the following we adopt as our central input value = 



1.43 GeV 1^ and use for the soft renormalization scale pLg = CFCis{lJ^s)'fnt- For rrit = 
174.3 GeV [H the latter gives /x, ^ 30 GeV. 

To evaluate the perturbative contribution we use the result of the previous section up 
to N^LO. In Figs. |^(a) and (b) the perturbative ground state energy in NLO, NNLO and 
N^LO is plotted as a function of the renormalization scale of the strong coupling constant 
for = and S" = 1, respectively. One can see, that the N^LO result shows a much 
weaker dependence on /i than the NNLO one. Moreover at the scale /i ~ 15 GeV for 
S = 1 and /i ~ 18 GeV for = 0, which are close to the physically motivated scale /i^, 
the N^LO correction vanishes and furthermore becomes independent of /x, i.e. the N^LO 
curve shows a local minimum. This suggests the convergence of the series for the ground 
state energy in the pole mass scheme, which is inconsistent with the renormalon picture. 
Transforming the result to the MS top quark mass leads to a worse convergence of the 
series for the resonance energy. However, the difference in the series behaviour due to the 
use of different mass parameters is less significant than in the case of the bottom quark 
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Figure 2: Ground state energy Ef'^' of the ti bound state in the zero- width approximation 
as a function of the renormahzation scale fi for 5* = (a) and = 1 (b). The short-dashed, 
long-dashed and solid line corresponds to the NLO, NNLO and N^LO approximations. 

as as for the top system is much smaller (see, e.g., [0). Note that the result shows very 
weak spin dependence. 

Collecting all the contributions we obtain a universal relation between the resonance 
energy in tt threshold production in e~^e~ annihilation or 77 collisions and the top quark 
pole mass 

mt - 174.3 GeV , , , 

Ere. = 1.9833 + 0.007 ± 0.0009 x , 21 

174.3 GeV 

where the weak nonlinear dependence of Eq. (|19D on mt is taken into account in the 
second term on the right-hand side of Eq. (^Tj). The central value is computed for mt = 
174.3 GeV. The sources of the theoretical uncertainty in Eq. (|2T|) are listed in Tab. |] 
in terms of the resonance energy. For the estimates we assume a ±100% error in the 
Fade approximation of the coefficient 03 and vary the normalization scale in the stability 
interval which is roughly given by 0.4 fig < ^ < fx^. Due to the very nice behaviour of the 
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perturbative expansion for the ground state energy we do not expect large higher order 
corrections to our result. 



4 Conclusions 

In this letter we have completed the result of [0 for the N'^LO corrections to the ground 
state energy of perturbative heavy quarkonium by computing the contributions propor- 
tional to the function, Eq. (|l^). We have found that the N^LO corrections are dominated 
neither by logarithmically enhanced a^ln(as) nor by the renormalon induced f^^a^ terms 
(cf. Eq. (|T^) ) and thus the full calculation of the correction is crucial for quantitative 
analysis. 

The result has been applied to determine the bottom quark mass from the mass of 
the T(IS') resonance. The perturbative series for the extracted bottom quark pole mass 
diverges and should be removed from the analysis in favour of the short- distance mass 
mbijnb) to get a reliable prediction. Our result (|^) for mbijrTi,) is in reasonable agreement 
with the NNLO estimates obtained in the framework of nonrelativistic T sum rules and 
the approach based on low moments, however, with a slightly higher central value. 



The explicit formula ([21| ) relating the top quark mass to the resonance energy in ti 
threshold production both in e~^e~ annihilation and 77 collisions has been derived. In 
contrast to the bottom quark case the perturbative expansion for the ground state energy 
of a (hypothetical) top-antitop bound state converges very well in the pole mass scheme. 
As a consequence the top quark pole mass can be extracted from the future experimental 
data on ti threshold production with the theoretical uncertainty of ±80 MeV. 
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